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ABSTRACT 

The  characteristics  of  general  laminar  flows  around  the  edge  of  a 
semi-infinite  flat  plate  are  investigated  on  the  basis  of  the  complete 
Navier-Stokes  equations.  The  solution  of  this  singular  and  nonlinear 
problem  is  obtained  in  an  exact  and  linear  manner  in  the  form  of  a 
double  series  around  the  edge.  It  permits  a  rigorous  discussion  of 
the  properties  of  flows  around  an  edge,  which  are  of  importance  in  the 
practical  applications  as  well  as  in  the  theory  of  viscous  and  non- 
viscous  fluid  motions.  In  particular,  it  becomes  possible  to  check 
the  reliability  of  the  Kutta  condition  which  represents  a  fundamental 
proposition  in  ideal  flow  theory. 


NWL  REPORT  NO.  1887 


FOREWORD 

This  work  was  sponsored  by  the  Naval  Weapons  Laboratory  Founda¬ 
tional  Research  Program.  It  was  performed  in  the  Computation  and 
Analysis  Laboratory  as  part  of  the  Meteorological  and  Oceanographical 
project  (R3 60FR103 / 2 10 1 /RO 1 10 1001 ) . 

The  authors  wish  to  extend  their  thanks  to  Dr.  C.  J.  Cohen  for 
many  critical  and  stimulating  discussions. 

The  date  of  completion  was  4  October  1963. 

Approved  for  Release: 

/s/  R.  H.  LYDDANE 

Technical  Director 


NWL  REPORT  NO.  1887 


1.  Introduction 

The  flows  around  sharp  edges  can  be  described  in  a  realistic  way 
only  if  friction  effects  are  included.  For  the  mathematical  treatment 
of  this  problem  the  complete  Navier-Stokes  equations  must  be  used. 
Early  attempts  to  solve  the  problem  on  the  basis  of  ideal  flow  theory 
by  the  introduction  of  discrete  vortices  or  continuous  vortex  sheets, 
are  in  principle  of  hypothetical  nature  and  must  be  confirmed,  a  pos¬ 
teriori,  by  the  real  flow  behavior.  It  turned  out  that  the  boundaries 
of  two-phase  flows  (for  instance,  cavities)  could  be  described  quite 
well,  whereas  the  approach  for  one-phase  flows  (for  instance,  wakes) 
is  not  satisfactory  (see  [2]). 

In  a  recent  paper  [7]  a  method  has  been  presented  to  solve  the 
edge  problem  with  the  Navier-Stokes  equations  for  the  special  case 
of  symmetric  flows  past  a  semi-infinite  plate.  The  same  procedure 
will  be  applied  in  this  paper  to  discuss  the  general  nonsymmetric 
flows  in  the  vicinity  of  the  edge.  Though  the  discussion  is  essen¬ 
tially  qualitative,  valuable  information  about  the  laminar  flow 
characteristics  are  obtained.  They  reveal  not  only  patterns  of  flows 
towards  or  from  an  oblique  plate  and  explain  phenomena  observed  in 
photographs,  but  they  permit  also  a  reexamination  of  propositions  in 
ideal  flow  theory.  For  example,  it  can  be  shown  that  separation  at 
the  edge  occurs,  which  contradicts  the  well-known  Kutta  condition. 
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2 .  Reduction  of  the  Equat ions  of  Motion 

For  a  two-dimensional  steady-state  flow  the  Navier -Stokes  equa¬ 
tions  and  the  equation  of  continuity  are  in  a  Cartesian  coordinate 
system  (x,  y) 


uux  +  vuy  ”  Px  +  u^uxx  +  ^y) 

UVX  +  Wy  =  -  i  Py  +  +  Vyy) 

Ux  +  Vy  =  0 


(1) 

(2) 

(3) 


where  (u,  v)  is  the  velocity,  p  the  pressure, 
and  u  the  constant  kinematic  viscosity  of  the 
Considering  the  flow  past  a  semi-infinite 
located  along  the  positive  x-axis,  one  has  to 
condition  at  the  plate 


p  the  constant  density, 
fluid, 

flat  plate,  which  is 
satisfy  the  nonslip 


(4) 


According  to  [7]  it  is  useful  to  remove  the  singularity  at  the 
origin  by  the  introduction  of  parabolic  coordinates,  which  are  made 
dimensionless  by  a  constant  length  a  .  This  may  be  chosen  in  such  a 
way  that  its  value  becomes  characteristic  for  the  flows  under 
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consideration  (see  [7]).  The  transformation  is  • 
x  =  o*(S2  -  Tf) 
y  =  2o2§ti  . 


(5) 

(6) 


The  velocity  and  the  pressure  are  also  replaced  by  dimensionless 
quantities 

U  =  —  u  ,  V  =  —  v  ,  P  =  -  2  —  .  (7) 

V  v  u3  p 

One  obtains  for  the  equations  (1)  through  (3) 

2(§u  +  T]V)U5  +  2(§V  -  Tp)^  =  -  (8) 

2(§U  +  71V)V§  +  2(§V  -  T]U)V^  =  T]P?  +  ^  (9) 

§(U§  +  V^)  +  T](V§  -  U^)  =0  (10) 

with  the  nonslip  condition 


T|  =  0  :  U  =  0 ,  V  =  0 


(ID 


at  the  plate.  The  analytic  properties  of  the  solutions  permit 
Taylor  expansions  of  the  integrals  around  §  =  0: 

u=  S  u„ CO) 5°  ,  v-  S  v„ (Tl) ?”  ,p  =  S  Pn(TD§°  (12) 

n=o  n=0  n— 0 
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With  these  series  the  partial  differential  equations  (8)  through  (10) 
are  reduced  to  the  infinite  set  of  ordinary  differential  equations 

<  +  (n  +  1)  (n  +  2)Un+s  -  T|Pn'  +  nP„  = 


n+1  n 

2  Z  iUjCUn-j  +  TVn-1+1)  +  2  Z  -  T|Ua- ± )  (13) 


i=0 

i=0 

V"  +  (n  +  1) (n  +  2)Vt+a  + 

(n  +  1)T1P»+1  +  P'-i  = 

n+l 

n  . 

2  S  iV,(Un-,  +  T|Vn_1+z  )  + 

2  r  V/CV.-,-!  -  TlUn_t) 

(14) 

1  =  0 

i=0 

nU„  +  (n  +  1)T|Vb+x  -  T|UB  + 

o 

II 

H 

*»  rt 
> 

(15) 

for  n  =  0,  1,  2  ...  .  The  primes  denote  the  derivatives  with 

respect  to  T|.  The  nonslip  condition  (11)  requires 

T|  =  0:  Un  =  0  ,  Vn  =  0  .  (16) 

For  the  sake  of  simplicity  in  the  forthcoming  derivations  the  first 
19  equations  of  the  system  (13)  through  (15)  are  listed  in  Appendix  A. 

The  coefficient  functions  Un ,  Vn ,  and  Pn  are  assumed  in  the  form 
of  series  expansions 

Un  =  S  U^lf  ,  V.  =  Z  Vn,Tf  ,  P»  =  T  P„.Tf  (17) 

m=l  1=2  B=0 

in  which  the  nonslip  condition  (16)  and  the  condition 

VBi  =  0  for  all  n  =  0,  1,  2  ...  (18) 
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have  been  considered  already.  The  latter  one  follows  immediately  from 
the  reduced  continuity  equations  (15).  As  was  pointed  out  in  [7]  the 
existence  of  an  integral  for  the  set  (13)  through  (15)  is  not  certain 
since  the  differential  equations  are  singular  at  the  origin  T)  =  0. 

The  existence  proof  is  achieved  in  the  same  way  as  in  [7]  by  substitut¬ 
ing  the  series  (17)  into  the  differential  equations  (13),  (14),  and 
(15).  These  reduce  to  an  infinite  set  of  algebraic  equations,  which 
consists  of  an  infinite  number  of  finite  subsets,  all  of  which  are 
linear  in  the  unknown  coefficients.  The  number  of  free  parameters 
for  each  linearly  dependent  subset  is  now  two  instead  of  one  parameter 
found  in  the  special  case  of  symmetric  flows.  They  are  designated 
by  An  and  Bn  in  such  a  way  that  all  Bn  vanish  for  symmetric  flows . 
Therefore,  the  family  of  all  symmetric  flows  may  be  defined  by  the 
condition 

Bn  =  0  (n  =  2,  3,  4,  ...)  .  (19) 

The  subsets  of  algebraic  equations  are  found  by  collecting  the 
algebraic  equations  of  the  same  order,  where  the  order  is  defined 
by  the  sum  of  the  exponents  of  7)  and  §  .  The  first  subset  yields 
vanishing  coefficients  because  of  the  conditions  (16)  and  (18). 

The  second  subset  follows  from 

(A  -  DTi1  :  Uos  =  0  (20) 


5 


NWL  REPORT  NO.  1887 


(A  - 

2)T?  : 

Un  -  Utl  =  0 

(21) 

(A  - 

9)  T|°  : 

O 

II 

TO 

O 

> 

(22) 

This  system  has  'the  solution 

U03  -  0  ,  Un  =  B2  ,  V02  =  0  (23) 

where  Ba  is  the  free  parameter.  The  pressure  is  normalized  by  assum¬ 
ing  a  zero  pressure  at  the  edge  of  the  plate: 

Poo  =  0  (24) 

The  next  subsets  up  to  the  seventh  order  are  recorded  in  Appendix  B. 
Their  two-parameter  solutions  are: 

Third  order  subset  (with  A3  and  B3  as  parameters) : 


U03  = 

A3 

V03 

=  b3 

Pol  “  “  6A3 

Ul3  = 

3B3 

Via 

=  3A3 

P10  =  -  6B3 

(25) 

U8i  = 

-  6A3 

Fourth  order  subset  (with  A4  and  B4  as  parameters) : 


U04  = 

0 

0 

11 

* 

0 

> 

P02  - 

A4 

Ul3  = 

b4 

to 

II 

0 

Pli  =  - 

2B4 

(26) 

11 

TO 

TO 

!=> 

a4 

v22  =  B4 

P20  =  - 

a4 

U3i  -  -  B4 
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Fifth  order  subset  (with  Ag  and  B5  as  parameters) : 


a 

o 

01 

II 

A5 

11 

10 

0 

> 

B5 

P03  =  -  10^5 

u14  - 

5B5 

Vl4  = 

5Ab 

Pis  -  -  30Be 

U23  =  “ 

25Ae 

V23  = 

53b 

Psi  =  30A 

U32  =  " 

15Bb 

V32  = 

-  15A5 

P30  =  10bb 

U41  = 

ioa5 

Sixth  order 

subset 

(with  Ag  and 

B6  as 

parameters) : 

U06  = 

0 

V06  = 

0 

P04  =  3A« 

& 

OT 

II 

b6 

Vl5  = 

0 

Pi  3  =  '  8Be 

U24  = 

6Aq 

V24  = 

4B6 

Pss  =  '  18As 

U33  =  ' 

6Be 

V33  = 

4Ag 

P31  =  8B6 

1 

II 

ct 

+ 

6Ag 

V42  = 

-  4B6 

P40  =  3Ag 

Ubi  = 

2B6 
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Seventh  order  subset  (with  A?  and  B7  as  parameters): 


<=: 

0 

11 

A? 

Vo  7  = 

B7 

Uie  = 

73,-1 

b2a3 

Vie  = 

a7 

I 

II 

ID 

n 

8A,  -  1 

B2B3 

VaB  = 

14B7  - 

1  T3  A 

2  B2A3 

U34  =  - 

70B,  +  | 

b8a3 

V34  =  - 

10A7  - 

^  B2B3 

U43  = 

15A,  +jf- 

b2b3 

V43  =  - 

35B7  + 

B2A3 

U62  = 

35B7  - 

b2a3 

v52  = 

5A7  + 

— bsb3 

14 


U6i  -  -  2A7  -  B2  B3 


Po5  *  -  2A7  -  ~  B2B3 

Pj  4  =  *  70B7  +  B2A3 

Ps3  s  20A7  +  ”  B2B3  (29) 

B32  =  140 B7  -  4B2A3 

P41  =  -  10A7  “  ~  B2B3 

Pbo  =  “  1^b7  +  B2A3 

With  the  equations  (12) ,  (17)  and  the  coefficients  (23)  through 
(29)  one  obtains  for  the  velocity  components  U,  V  and  the  pressure  P 
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the  double  series 
U  *  B27)5 

-  6A3T)?8.+  3Bslf|  +  A3lf 


-  B4-TI53  +  A^l3  +  B4lf  | 

+  lOAeTl4  -  XSBsTl®?3  -  25A6T|3Ss  +  5BsTf  5  +  As  if 
+  2 Belli6  -  SAeH3!4  -  6Belf|3  +  6Aelf|3  +  BBlf  5 

-  (2A7  +  JJB2B3)T1|8  +  (35B7  -  B2A3)Tp|B  +  (15AV 

-  (70B,  -  §B2A3)Tf|3  -  (8A7  +  7B2B3)lf|3  +  (7B, 

1  „7 

+  7^7  T| 


+  B2B3)Tf  |4 

-  3B2A3)f  | 

(30) 


V  =  3A3Tp|  +  B3Tf 

-  B4Tf|3 

-  15A,Tf§3  +  5BETf?8  +  5Aslf  |  +  Be  if 

-  4B6Tf|4  +  4Aelf|3  +  4B8  if  I3 

+  (5 A,  +^B2B3)Tf|g  -  (35B,  -  B2A3)lf|4  -  (10A7  +  )b2B3)1T4|3 
+  (14B,  -  7B2A3)Tf?a  +  A,  If  I  +  B,lf 


P  =  -  6B3?  -  6A3T] 

-  A*!3  -  2B4HI  +  A4lf 

+  10Bg|3  +  30 A5 1113  -  SOBglfS  -  10Aslf 

+  3  As  I4  +  8BS  1113  -  18Ae  If  I3  -  8Belf|  +  3Ag  if 

-  (14B7  -  •JB2A3)|6  -  (10A,  +  ^B2B3)T||4  +  (140 B7  -  4B2A3)lf|3 

+  (20A7  +  |b2B3)H3|3  -  (70B,  -  BaA3)lf  |  -  (2A7  -  |^B2B3)lf  (32) 


9 


NWL  REPORT  NO.  1887 


The  description  of  the  flow  behavior  by  the  first  terms  of  the 
series  expansions  (30) ,  (31)  ,  and  (32)  is  valid  in  the  neighborhood 
of  the  edge.  It  is  significant  that  in  these  double  series  nonlinear 
inertial  terms  do  not  appear  before  the  seventh  order.  Therefore, 
the  flow  region  near  the  edge  behaves  like  a  flow  field  of  Stokes' 
type.  This  is  actually  the  key  to  the  understanding  of  the  detailed 
description  of  the  flow  properties  without  knowing  the  asymptotic 
flow  far  away  from  the  edge.  As  was  pointed  out  in  [6]  and  [7], 
every  flow  forms  a  "s low-mot  ion  layer"  at  a  rigid  body,  for  which  the 
inertial  forces  may  be  neglected.  However,  the  inertial  forces  become 
dominant  outside  the  s low-mot ion  layer  and  determine  the  asymptotic 
ideal  flows.  For  this  reason,  a  slow-motion  flow  is  in  principle  not 
possible,  if  it  is  required  to  have  asymptotic  ideal  flow  properties. 

In  this  connection  a  special  case  of  the  general  solutions  (30) , 
(31),  and  (32)  may  be  anticipated,  which  belongs  to  the  family  of  the 
Couette  flows  and  which  is  characterized  by  the  property  that  the 
plate  does  not  disturb  the  surrounding  flow.  Assuming  all  coefficients 
to  be  zero  except  Bs  and  A4  one  obtains 

U  =  Ban?  +  A.,  If  Is  (33) 

V  =  0  (34) 

P  =  A* (If  -  Is)  (35) 
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or 


u 


ifk 

7W 


y  + 


v  =  o 


(36) 

(37) 


p  =  A4 

p  2a4  “7 


(38) 


In  contrast  to  the  symmetric  Couette  flows,  which  are  defined  by 
B2  =  0,  the  nonsymmetric  flows  (36),  (37),  and  (38)  cause  a  shear 
stress  at  the  plate  of  the  magnitude 

Upu  (x,  0)  =  p2p  ^4  ,  (39) 

y  2a 

which  has  a  positive  sign  at  the  upper  side  of  the  plate  and  a  nega¬ 
tive  sign  at  the  lower  side,  so  that  they  balance  each  other  and 
yield  no  resistance. 

It  might  be  mentioned  that  the  existence  of  the  slow  motion  layer 
at  a  rigid  body  is  one  reason r  why  Prandtl’s  boundary  layer  theory 
cannot  describe  the  "thin  friction  layer"  at  the  body  in  an  exact 
manner  since  the  friction  terms  in  the  Navier-Stokes  equations  are 
truncated  (see  [6],  [7]). 

3 .  Discussion  of  Flow  Properties 

Each  combination  of  the  free  parameters  An  and  Bn  represents  an 
exact  solution  of  the  Navier-Stokes  equations.  From  the  infinite 
number  of  possibilities,  which  reflects  the  generality  of  the 
Navier-Stokes  equations,  only  simple  cases  of  practical  interest  may 
be  pointed  out.  These  selected  examples  may  be  divided  into  attached 
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and  separated  flows  around  the  edge. 

The  distinction  between  these  two  classes  follows  from  the  study  of 
the  streamlines  at  the  edge.  If  no  streamlines  originate  at  the  edge, 
the  flow  may  be  called  attached;  otherwise  the  flow  is  separated.  The 
latter  class  includes  also  the  trivial  cases  such  as  the  parallel  flows 
toward  or  away  from  the  plate  (Figures  la  and  2a) . 

The  streamlines  at  the  edge  are  described  by  the  differential  equa¬ 
tion 

dy  7)d§  +  §d7]  v  V 


dx  §d§  -  TJdT|  u  U 


(40) 


or 


£5  -  BaH5  +  A3H(95a  -  Tf)  +  ... 

d§  "  B25a  +  As?(4Tf  -  6?a)  +  ... 


(41) 


In  addition,  it  is  useful  to  know  the  shear  stress  at  the  plate,  which 
follows  from  equation  (30) 


vpuv(x,  0)  =  UflCS,  o)  =  -^(Bs  -  6A3|  -  B4?2  +  10As?3  +  ...)  (42) 

y  2(7  §  1  2a 

The  discussion  of  the  flow  properties  may  be  guided  by  envisioning  the 
following  experiment:  Let  the  motion  be  symmetrical  around  the  semi¬ 

infinite  plate,  where  the  flow  toward  the  plate  is  called  stagnation 
flow  past  the  (leading)  edge  and  the  flow  away  from  the  plate  is 
designated  as  wake  flow  around  the  (trailing)  edge  (Figures  la  and 
2a) .  Keeping  the  Reynolds  number  constant  one  may  change  the  flow 
direction  slowly  up  to  a  right  angle  of  attack.  One  observes  then 


u8p 
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the  following  flow  behavior:  Under  a  sufficiently  small  angle  of  attack 
a  streamline  will  always  end  or  begin  at  the  edge  exactly  as  in  the 
symmetric  case  (Figures  lb  and  2b) .  For  larger  angles  of  attack  this 
streamline  vanishes  from  the  edge,  and  the  fluid  particles  pass 
smoothly  around  the  leading  and  trailing  edges  (attached  flows  in 
Figures  lc  and  d,  2c  and  d) .  The  flow  remains  attached  as  long  as  the 
Reynolds  number  is  sufficiently  small.  Otherwise,  from  a  certain  angle 
of  attack  on  separation  occurs  (Figures  3  and  4). 

From  this  description  the  following  types  of  flows  are  distinguished 

a.  Symmetric  flows 

They  are  discussed  already  in  [7]  and  are  defined  by  the  condition 

Aa  =  0,  A4A5  >  0  ,  Bn  =  0  (n  =  2,  3,  ...)  (43) 

The  flows  around  leading  and  trailing  edges  are  determined  by  the 
positive  and  negative  signs  of  the  coefficient  A5 .  Condition  (43) 
includes  only  symmetric  motions  of  the  trivially  separated  flow  type. 

b.  Quasi-symmetric  flows 

In  contrast  to  ideal  flow  theory,  where  an  infinitesimal  small 
angle  of  attack  causes  the  disappearance  of  the  dividing  streamline 
from  the  edge,  the  real  flow  can  stand  a  slight  nonsymmetry  without 
losing  the  properties  of  symmetric  flows  in  the  leading  terms.  The 
conditions  are  the  same  as  in  (43)  except  that  a  nonsymmetry  arises 
from  the  seventh  order  on  (Figures  lb  and  2b) : 

A3  =  0,  A4A5  >0.  B„  =  0  for  n  <  7,  B,  /  0  .  (44) 
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The  fact  that  a  streamline  at  the  edge  exists  is  anticipated  and  will 
be  shown  farther  below. 

Because  of  the  same  behavior  of  symmetric  and  quasi-symmetric  flows 
the  theorems  established  in  [7]  for  symmetric  flows  are  valid  for  both 
classes.  Therefore,  one  obtains  for  the  shear  stress  and  the  pressure 
the  following  extended  theorems : 

THEORfM  1:  In  symmetric  and  quasi-symmetric  laminar  flows  past  a  semi¬ 
infinite  flat  plate  the  shear  stress  at  the  surface  is  zero  at  the  edge 
and  increases  or  decreases  as  the  three-halves  power  of  the  distance 
from  the  edge. 

THEORfM  2:  In  symmetric  and  quasi- symmetric  laminar  flows  past  a  semi¬ 
infinite  flat  plate  the  pressure  varies  at  the  edge  in  a  strictly  linear 
manner. 

c.  Almost  skew- symmetric  flows 

For  larger  angles  of  attack  the  streamline,  which  originates  at  the 
edge,  vanishes,  and  the  flow  direction  at  the  upper  and  lower  side  of 
the  plate  is  opposite.  The  coefficient  B2  is  no  longer  zero.  However, 
one  observes  still  a  symmetric  property  which  is  caused  by  the  coef¬ 
ficient  B4  in  equation  (31).  This  type  of  motion  may,  therefore,  be 
called  "almost  skew-symmetric'’  flow  and  is  described  by  the  condition 

A3  =  0,  A5  >  0,  B2B3  >  0,  B3B4  >  0  (45) 

for  the  leading  edge  and 

A3  =  0,  A5<  0,  B2B3  >  0  ,  B3B4  <  0  (46) 

for  the  trailing  edge  (Figures  lc  and  2c). 
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The  non-zero  coefficient  Bs  results  in  a  change  of  the  flow  charac¬ 
teristics  which  become  similar  to  the  type  of  Couette  flows  (33), 

Indeed,  the  lowest  order  of  the  velocity  component  u  is  linear  in  the 
coordinate  y.  The  shear  stress  at  the  plate  is  constant  in  the  very 
neighborhood  of  the  edge.  This  result  differs  from  that  for  symmetric 
and  quasi-symmetric  flows  formulated  in  Theorem  1. 

The  pressure  drop  at  the  plate  is  also  different  from  that  of  the 
symmetric  and  quasi-symmetric  flows  (Theorem  2),  As  the  coefficient 
B3  is  not  zero  the  pressure  becomes 

P  =  P  ±6  h.  Jk  -  ^  x  ±  ...  \  (x  ^  0,  y  =  0)  (47) 

2<j  \  cr  o*  / 

Therefore,  the  pressure  decreases  near  the  edge  in  flow  direction  with 
the  square  root  of  the  distance  from  the  edge.  This  fact  is  surpris¬ 
ing  from  the  standpoint  of  ideal  flow. theory,  where  an  extremum  of  the 
body  contour  effects  an  extremal  pressure,  A  comparison  of  the  pres¬ 
sure  distribution  around  a  circular  cylinder  for  ideal  and  real  flows 
shows  on  the  other  hand  [5]  that  tne  expected  extremal  behavior  is  not 
necessarily  true  for  real  flows, 
d.  Distorted  nonsymmetric  flows 

A  further  increase  of  the  angle  of  attack  leads  to  a  type  of  flows 
which  has  already  the  tendency  to  separate  (Figures  Id  and  2d) .  This 
is  due  to  the  fact  that  the  coefficient  A3  is  no  longer  zero.  One  can 
easily  see  that  the  term  SAgTpS  in  equation  (31)  for  the  velocity  com¬ 
ponent  V  distorts  the  symmetry. 
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The  flows  of  this  type  are  characterized  by 

A3  £0  ,  Ba B3  >  0  (48) 

for  the  leading  and  trailing  edges,  respectively. 

The  statements  on  shear  stress  and  pressure  drop  at  the  plate  for 
almost  skew- symmetric  flows  are  also  valid  for  distorted  nonsymmetric 
flows.  Therefore  general  nonsymmetric  flows  have  the  following  prop¬ 
erties: 

THEOREM  3:  In  almost  skew-symmetric  and  distorted  nonsymmetric  laminar 
flows  around  the  edge  of  a  semi-infinite  flat  plate  the  shear  stress  at 
the  surface  is  constant  near  the  edge.  The  flow  behaves  like  a  Couette 
flow. 

THEOREM  4:  In  almost  skew-symmetric  and  distorted  nonsymmetric  laminar 
flows  around  the  edge  of  a  semi- infinite  flat  plate  the  pressure  at  the 
surface  decreases  in  flow  direction  as  the  square  root  of  the  distance 
from  the  edge, 
e.  Separated  flows 

Distorted  nonsymmetric  flows  around  the  edge  can  lead  to  separation 
from  the  edge.  Separated  flows  of  this  type  are  characterized  by  the 
existence  of  nontrivial  solutions  of  the  differential  equation  (41). 

A  solution  is  assumed  by  the  series  expansion 

§  =  all  +  eif  +  ylf+  ...  ,  (49) 

which  reduces  the  differential  equation  (41)  to  an  infinite  set  of 
algebraic  equations.  The  first  four  of  these  equations  are: 
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Tf  :  Baa*  =  0  (50) 

T)3  :  I5kac?  -  5B3q^  -  5A3Q'  -  B3  =  0  (51) 

Tf  :  3BS32  +  12B3q,3  +  6(1  -  9o2)A33  +  3(1  -  q?)q? B4 

+  2A4CV3  =  0  (52) 

“if  :  35A6q,b  -  35Bsa4  -  7(10A5  +  23B4)o3 

+  7 (2Bb  -  9vA3  +  3A4)c^  +  7(2yB3  -  9A332  +  Ag 
+  B43)ck'  +  7A3y  +  7B3(32  +  7Ba3y  +  Be  =  0  (53) 


There  follows  immediately  from  equation  (50)  that  a  non- vanishing 
coefficient  a  requires  the  coefficient  B2  to  be  zero.  In  the  physical 
(x,  y)  -  plane  the  condition  a  ^  0  means  that  the  plate  is  not  tan¬ 
gential  to  this  streamline.  In  this  case  motions  which  display  Couette 
flow  properties  (Ba  ^  0)  are  not  possible. 

It  is  easy  to  verify  by  means  of  the  velocity  components  U  and  V 
that  under  the  assumption  A3B3  /  0  only  the  combination 

Aa  >0  ,  >  °  (54) 

is  of  practical  interest.  This  class  of  flows  is  identified  as 
separated  flows  around  the  leading  edge  of  a  plate  (Figure  3). 

The  shear  stress  at  the  surface  is  now  according  to  equation  (42) 

J3  — 

Upuy(x,  0)  =  -  3  A3  *Jk  -  ...  (55) 

on  both  sides  of  the  plate. 


17 


NWL  REPORT  NO.  1887 


THEOREM  5:  In  separated  laminar  flows  around  the  leading  edge  of  a 
semi- infinite  flat  plate,  where  the  plate  is  not  tangential  to  the 
separating  streamline,  the  shear  stress  at  the  surface  is  zero  at 
the  edge  and  increases  as  the  square  root  of  the  distance  from  the 
edge. 

The  pressure  drop  at  the  plate  is  given  by  equation  (47)  and  is 
the  same  as  for  nonsymmetric  attached  flows.  However,  the  flow 
direction  is  in  the  present  case  the  same  on  both  sides  of  the  plate. 

The  coefficient  a  is  determined  by  the  cubic  equation  (51)  which 
shows  immediately  the  existence  of  a  real  root  and  which  may  be 
written  in  the  form 

3A3(o's  -  5)  =  ^(a2  +  f)  .  (56) 

Considering  the  restriction  (54)  one  realizes  the  condition 

(57) 

which  results  in  limiting  angles  2tt/3  and  -  2tt/3  for  positive  and 
negative  coefficients  B3,  respectively.  A  stronger  condition  for  a 
can  be  obtained  by  excluding  the  appearance  of  two  zero  points  of  the 
velocity  component  U  near  the  edge.  Then,  the  angle  between  plate 
and  streamline  at  the  edge  must  be  equal  or  smaller  than  —  . 

The  described  type  of  separation  at  the  leading  edge  can  be 
observed  on  photographs  [8,  page  683],  However,  there  results  an 
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important  consequence  for  the  ideal  flow  theory  which  will  be  explained 
briefly. 

In  ideal  flow  theory  vorticity  is  introduced  either  in  form  of  dis¬ 
crete  vortices  or  in  form  of  continuous  vortex  sheets  in  order  to 
remove  the  singularity  in  the  velocity  at  the  edge  of  a  body.  The 
strengths  of  these  vortices  or  the  vortex  density  of  the  sheets  are 
determined  by  the  well-known  Kutta  condition  which  requires  a  finite 
velocity  at  the  edge  (see,  for  instance  [l]).  This  condition  implies 
a  tangential  streamline  at  the  edge.  Also  Kirchhoff's  discontinuity 
lines  (free  streamlines)  are  tangential  at  the  edge  (see,  for 
instance  [2]).  These  propositions  are  obviously  in  contradiction 
to  the  real  flow  behavior,  at  least  for  one-phase  flows.*  It  is 
interesting  that  Rott  [4]  already  questioned  the  validity  of  the  Kutta 
condition  for  a  single  vortex  model.  He  found  an  angle  of  separation 
of  exactly  ^  on  the  basis  of  ideal  flow  theory  by  abandoning  the  Kutta 
condition. 


Equation  (50)  offers  besides  a  /  0,  B2  =0  also  the  alternative 
Of  =  0,  B2  /  0.  The  equations  (51),  (52),  and  (53)  yield 


B3  =  0 

(58) 

Pi  =  o,  p2  =  -  2  ti 

(59) 

b2 

Yl  -  -  If-  ,  ya  =  +  h. 

(60) 

7A3  7A3 

respectively. 

*  In  cavities  the  assumption  of  tangential  discontinuity  lines  is 
realistic  since  the  tangential  velocity  component  dominates  at  the  edge 
in  the  presence  of  a  two-phase  boundary. 
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For  non-vanishing  coefficients  A3  and  Bs.  two  streamlines  originate 

at  the  edge  under  180°  with  curvatures  which  can  be  expressed  by  the 

coefficients  =  0  and  $2  =  -  2  —  .  An  inspection  of  the  flow  patterns 

B2 

by  means  of  the  velocity  components  U  and  V  reveals  that  only  the  class 

Aa  <  0  ,  B2BS  >  0  (61) 

is  of  interest  for  the  present  discussion  (Figure  4).  This  class  may 
be  called  "separated  flows  around  the  trailing  edge  of  a  plate." 

The  shear  stress  at  the  surface  is  now  (because  of  B2  ^  0)  of 
Couette  type.  The  pressure  drop  is  the  same  as  for  the  separated  flows 
around  the  leading  edge. 

The  separated  flow  described  near  the  edge  (Figure  4)  can  be  part 
of  a  flow  field  which  is  produced  by  a  source  located  at  the  upper 
right  side  of  the  semi-infinite  plate  and  which  is  sketched  in  Figure  5. 
This  indicates  the  occurrence  of  "secondary"  vortices  near  the  edge. 

Rott  [4]  has  pointed  to  photographs  which  clearly  show  the  existence 
of  secondary  vortices  despite  their  tiny  extension  (see,  for  instance, 
the  flow  pattern  near  a  wedge  by  Prandtl  in  [3]). 

Here  again  the  propositions  of  ideal  flow  theory  must  be  reexamined. 
One  can  immediately  see  that  a  single  vortex  model  or  a  vortex  sheet, 
which  is  either  continuous  or  composed  of  discrete  vortices,  are  not 
an  adequate  means  to  describe  this  type  of  separation. 
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It  remains  to  check  the  existence  of  a  streamline  at  the  edge  which 
justifies  the  introduction  of  the  quasi -symmetric  flows.  Adopting  the 
conditions  A3  =  Bs  =  or  =  0  from  the  symmetric  trivially  separated  flows 


one  finds  with  the  equations  (51)  through  (53): 

B3  =  0  (65) 

Bsg3  =  0  (66) 

B6  •  0  (67) 

and  from  the  equations  of  the  next  higher  orders 

A53  =0  (68) 

3A43s(6  +  3)  +  9A5y  +  B7  =  0  (69) 


Since  AB  ^  0  there  follows  3  ■  o  and 


(70) 


The  coefficient  therefore,  causes  the  slight  nonsymmetry  of  the 
quasi-symmetric  flow  as  was  anticipated. 

It  may  be  mentioned  that  the  double  series  (30)  ,  (31)  ,  and  (32) 
for  the  velocity  components  U,  V  and  the  pressure  P  permit  the 
existence  of  other  classes  of  either  attached  or  separated  flows 
past  the  edge  of  a  plate.  For  instance,  the  lowest  order  terms  might 
vanish  and  higher  order  terms  become  of  leading  importance. 
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Also,  if 

A3  >  0  ,  B2B3  <  0  (71) 

the  differential  equation  for  the  streamlines  (41)  yields  a  solution  of 
the  form 

/  2  B3  J 

5  =  ±  Vs  bT  ^  +  •••  (72) 

which  is  not  analytic  at  the  origin  7]  =  0  as  is  assumed  in  equation  (49) . 
However,  a  simple  examination  of  the  flow  directions  reveals  that  these 
flows  are  of  only  remote  practical  interest  and  may,  therefore,  be  omitted 
in  this  discussion. 

Finally,  it  may  be  mentioned  that  an  introduction  of  polar  coor¬ 
dinates  around  the  edge  of  the  plate  would  have  led  to  exactly  the 
same  results.  This  method  will  be  presented  in  a  forthcoming  paper 
which  deals  with  the  more  general  wedge  and  corner  problems  that 
necessitate  the  introduction  of  polar  coordinates. 
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THE  LEADING  REDUCED  EQUATIONS  OF  MOTION 


a.  The  reduced  equations  of  continuity 

n  =  0  : 

Vi  -  U0'  =0 

n  =  1  : 

Ui  +  271  Va  -  1)  Ui  +  Vo  -  o 

n  =  2  : 

2U2  +  371  V3  -  71  Ug  +  Vi  =  0 

n  =  3  : 

3U3  +  471  V4  -  71  Ua  +  V2  =  0 

n  =  4  : 

4U4  +  511  Vs  -  71  U4  +  V3  =  0 

n  =  5  : 

5UB  +  6T]  V6  -  71  Ug  +  V4  =  0 

n  =  6  : 

6U6  4-  771  V7  -  1]  U6'  +  Vg  =  0 

b.  The  reduced  first  Navier-Stokes  equations: 
n  =  0  :  uj  +  2U8  -  T)  P0'  =  2T|  -  2TJ  U0U0' 

n  =  1  :  Ui  +  6U3  -  T)  P*  +  P:  = 

2U0Ui  +  4T)  V0U2  -  2T]  U0Ui  +  2V0U0' 

n  =  2  :  Ug  +  12U4  -  1]  Ps'  +  2P2  = 

2Uf  +  271  UiV2  +  4U0U2  +  271  U0'u2  +  67]  U3V0 
+  2Ui'v0  -  21]  Ui'Ux  +  2Uq2~  271  U2U0 
n  =  3  :  Ug  +  20U5  -  T|  P3  +  3P3  - 

eUxUg  +  271  Ui  V3  +  411  U2V2  +  6U0U3  +  411  U3U0' 

-  811  U4V0  +  2U0'vs  +  2UiUq  -  2J\  u/u2  +  2UgV0 

-  21]  UgUr  -271  U3U0 


(A  -  1) 
(A  -  2) 
(A  -  3) 
(A  -  4) 
(A  -  5) 
(A  -  6) 
(A  -  7) 

(A  -  8) 

(A  -  9) 

(A  -  10) 


(A  -  11) 
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n  =  4  :  U4  +  30US  -  T)  P4'  +  4P4  = 

8UiU3  +  27]  U1U4  +  4U§  +  47]  U2V3  +  671  U3V2  +  8U0U4 
+  67]  U4Uq  +  107]  U5V0  +  2UqV3  +  2V2Ui  -  27]  U3u/ 

+  2UoUs  -  27]  UgUg  +  2V0U3  -  27]  U1U3/  -  27]  U0U4'  (A  -  12) 
n  =  5  :  Ug  +  42U7  -  7]  PB'  +  5PB  = 

2V4Uq  +  87]  U5U0  +  2V3Ui  -  27]  U4u/  +  2U4Ui 

+  27]  V5U0  +  2VsUa  -  27]  U3U2  +  4U3U2 

+  47]  V4U2  +  2UqU3  -  27]  U2U3  +  6U2U3  +  67]  V3U3 

+  2V0U4  -  27]  U1U4/  +  8U1,  U4  +  87]  V2U4  -  27]  U0Ub 

+  10UoUB  +  127|  V0U6  (A  -  13) 

c.  The  reduced  second  Navier-Stokes  equations: 

n  =  0  :  Vo  +  2VS  +  11  Pi  =  211  V0U0'  -  2T|  U0V0'  (A  -  14) 

n  =  1  :  Uq  +  6Va  +  2T)  Ps  +  P0'  = 

2U0'u0  +  41)  V0Vg  -  211  U0Uo  +  211  U0/2 

+  2V0'V0  -  211  VqUj  (A  -  15) 

n  =  2  :  Vs  +  12V4  +  31]  P3  +  P/  = 

2U0'U!  +  611  U0'Vg  +  4ii0Vg  -  611  V0V3  +  2U0V0' 

-  211  V0'Ug  +  2V0U"  -  2H  u£Uj  -  2T|  U0Vg'  (A  -  16) 

n  =  3  :  Vg  +  20 Vg  +  41)  P4  +  Pg  = 

2U0'Ug  +  81]  U0'V3  +  4Ux  V2  +  411  Vl  +  6U0V3 

-  811  V0V4  +  2V3Vg  -  211  V0'U3  +  2UoU0'  -  211  U^Ug 

+  2V0Vg'  -  211  UjVg'  -  211  U0V3  (A  -  17) 
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n  =  4  :  V'l  +  30 Vs  +  5T|  PB  +  P3'  + 

2UaU0'  +  XOT)  U0'v4  +  4U2V2  +  1011  V2V3 
+  6Ui  V3  +  8U0V4  +  1011  V0VS  +  2V3Vq 

-  2H  U4V0'  +  2V2U3  -  211  U3U"  +  2U0'V2 

-  211  U2V2'  +  2V3V3  -  211  Uj.V3'  -  211  U0V4  (A  -  18) 

n  =  5  :  Vg'  +  42V,  +  611  P6  +  P4  = 

2V4V3  -  2T|  UBV0'  +  2U4Uq  +  1211  VBU0' 

+  2V3U3  -  2T|  U4U"  +  4U3V2  +  41|  V2V4 

+  2V2V2  -  211  U3V2  +  6U2V3  +  611  Vi 

+  2U0'V3  -  211  U2V3  +  8Ui  V4  +  811  V2V4 

+  2V0V4  -  211  UXV4  +  10UoVB  -  211  U0VB  +  1211  V0V6  (A  -  19) 

The  equations  (A  -  8)  through  (A  -  19)  are  simplified  by 
considering  equation  (A  -  1). 
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THE  FIRST  SEVEN  SUBSETS  OF  ALGEBRAIC  EQUATIONS. 

First  subset: 

Vanishing  coefficients  because  of  conditions  (16)  and  (18) . 


Second  subset: 
From  equations 


(A  -  1),  I)1  : 

Uos  =  0 

(A  -  2) ,  Tf  : 

Un  -  Uu  =  0 

(A  -  9)  ,  Tf  : 

vos  =  o 

Third  subset: 

From  equations 

(A  -  2) ,  Tf  : 

3B3  -  Ula  =  0 

(A  -  3) ,  Tf  : 

U31  +  6A3  =  0 

/-s 

> 

00 

'-Z 

2Uaj  +  6A3  -  P0i  =  0 

(A  -  14), Tl1  : 

6»3  +  Pio  “  0 

(A  -  9) ,  T|°  : 

2Uis  +  P10  =  0 

(A  -  15), Tf  : 

*oi  +  2Vla  =  0 

Fourth  subset: 

(A  -  2),  Tf:  Ula  +  4B*  +  2Vaa  -  3U13  =  0 
(A  -  8),  Tf:  -  2Poa  +  2Uaa  +  12A*  -  0 

(A  -  14),  Tf:  PU  +  2Vsa  +  12B4  =  0 

(A  -3),  Tf:  A*  =  0 
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(A  -  9)  ,  Tf  : 

U3i  +  Uls  =  0 

(A  -  15),  T)1  : 

2PS0  +  2Pos  +  6Vj.  3  =  0 

(A  -  4)  ,  Tf  : 

U»l  +  v2s  =  0 

(A  -  10),  T|°  : 

2PS0  T"  2U2a  =  0 

(A  -  16),  T|°  : 

P11  +  2VSS  =  0 

Fifth  Subset 

(A  -  2)  ,  Tf  : 

5BS  +  2Vgs  -  3U14  =  0 

(A  -  8),  Tf  : 

-  3P03  +  2U23  +  20A6  =  0 

(A  -  14) ,  Tf  : 

Pis  +  2VS3  +  20  B5  =  0 

(A  -  3)  ,  Tf  : 

20AS  -  US3  +  3V32  =  0 

(A  -  9),  Tf  : 

-  Pi  a  +  6U33  +  120i4  =  0 

(A  -  15),  Tf  :  - 

2Pgi  +  3P03  +  6V3S  +  12V!  4  =  0 

(A  -  4),  Tf  : 

U3S  +  3Vaa  =  0 

(A  -  10),  Tf  : 

P21  +  12U4i  +  6Usa  =  0 

(A  -  16),  Tf  : 

3Pao  +  2Pu  +  6Vsa  =  0 

(A  -  5) ,  Tf  : 

3U4i  +  2vS2  =  0 

(A  -  11) ,  Tf  : 

3P30  +  2U32  =  0 

(A  -  17),  T|°  : 

Pai  +  2V32  =  0 

Sixth  Subset 

(A  -  1) ,  Tf  : 

V1S  -  6As  =  0 

(A  -  2)  ,  Tf  : 

6BS  +  2VS4  -  4Ui5  =  0 

(A  -  8) ,  Tf  : 

-  4P04  +  2U34  +  30Ae  =  0 

(A  -  14) ,  If  : 

Pi  3  +  2VS4  +  30B6  =  0 
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(A  -  3), 

Tf  : 

30AS  +  3V33  -  2U24 

(A  -  9) , 

Tf  : 

-  2Pi3  +  6U33  +  20Ui  5 

(A  -  15)  , 

7]3  : 

2P2S  +  4P04  +  6V33  +  20Vis 

(A  -  4)  , 

'If  : 

V34  +  V4g 

(A  -  10), 

: 

Ug4  +  U4g 

(A  -  16)  , 

r?  : 

P31  +  Pi  3  +  4V42  +  4Vg4 

(A  -  5)  , 

Tf  : 

4U42  +  3V33  -  2U4g 

(A  -  11)  , 

T]1  : 

2P3i  +  20US1  +  6U33 

(A  -  17), 

Tl1  : 

4P40  +  PPss  P  ^33 

(A  -  6), 

'll1  •* 

2Usl  +  V4g 

(A  -  12)  , 

Tl°  : 

2P40  +  U4g 

(A  -  18)  , 

T]°  : 

P3l  +  2V*g 

=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 
=  0 


Seventh  Subset 


(A  - 

1), 

Tf  : 

7A,  -  VlR 

= 

0 

(A  - 

2), 

Tf  : 

-  5U1S  +  7B,  +  2VgS 

= 

0 

(A  - 

8), 

Tf  : 

-  5P0S  +  2U25  +  42A7 

= 

2B3Ui  1 

(A  - 

14), 

Tf  : 

pl4  +  2VgS  +  42 B7 

= 

0 

(A  - 

3), 

Tf  : 

3V34  -  3UgS  +  42A, 

= 

0 

(A  - 

9), 

Tl4  : 

-  Pj  4  +  2U34  +  1014  s 

= 

0 

(A  - 

15), 

Tl4  : 

2Pg3  +  5P0S  +  6V34  +  30Vi6 

= 

-  6B3Uh 

(A  - 

4), 

Tf  = 

5VSS  +  4V43  -  U34 

= 

0 

(A  - 

10), 

Tf  : 

-  P23  +  I2U43  +  20UgS 

= 

2B3Ui 1  - 

2U:  2Ui  i 
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(A  - 

16), 

if  : 

3P32 

+  4P14  + 

12V43  + 

20V35 

=  6A3Uh 

-  4VlsUn 

(A  ■ 

5), 

Tf  : 

U43 

+  4V34  +  5Vbs 

=  0 

(A  ■ 

ID, 

if  : 

P32  + 

20Ub3  + 

12U34 

=  2UaiUu 

+  2VlsUn 

(A  - 

■  17), 

if  : 

4P4i 

+  3P23  + 

20VS2  + 

12V34 

=  0 

(A  - 

•  6), 

If  : 

Ub2 

+  v13 

=  0 

(A  - 

12), 

Tl1  : 

P41  +  10U6i 

+ 

2U43 

=  0 

(A  - 

18), 

Tl1  : 

31*60  2P32 

+ 

6V43 

=  0 

(A  ■ 

7), 

Tl1  : 

5U61 

+ 

2V6S 

=  0 

(A  - 

•  13), 

Tf  : 

0 

10 

PL, 

m 

+ 

2U6e 

=  0 

(A  - 

'  19), 

Tf  : 

P41 

+ 

2V6e 

=  0 
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*■ 


A3=0 ,  Ki>  0,  Ab>  i> 
Bn  =  0 


X 


a3=o,  a4>  0,  Ab>  0 
Bn=0  for  n  <  7 , 


B7>  0 


A3=0,  Ab  >  0, 
B2B3>  0,B3B4-'>  0 


A3>  0, 

x-  B2B3  >  0 


leading  edge 
angles  of 


:  Sketch  of  a  separated  flow  around  the  trailing 
edge  of  a  semi- infinite  flat  plate  with  a  source 
on  the  upper  right  side. 
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